The quantum adiabatic algorithm and scaling of gaps at first order quantum phase transitions 
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Motivated by the quantum adiabatic algorithm (QAA), we consider the scaling of the Hamiltonian gap at 
quantum first order transitions, generally expected to be exponentially small in the size of the system. However, 
we show that a quantum antiferromagnetic Ising chain in a staggered field can exhibit a first order transition 
with only an algebraically small gap. In addition, we construct a simple classical translationally invariant one- 
dimensional Hamiltonian containing nearest-neighbour interactions only, which exhibits an exponential gap at 
a thermodynamic quantum first-order transition of essentially topological origin. This establishes that (i) the 
QAA can be successful even across first order transitions but also that (ii) it can fail on exceedingly simple 
problems readily solved by inspection, or by classical annealing. 

PACS numbers: 03.67.Ac, 05.30.Rt, 75.10.Jm, 75.10.Pq, 75.30.Kz 



The quantum adiabatic algorithm (QAA) |[T]|2l promises to 
harness the power of quantum mechanics, in particular quan- 
tum tunneling through energy barriers, in order to solve hard 
optimization problems more efficiently. Like classical simu- 
lated annealing (CSA), part of its appeal is its general purpose 
"black box" nature. Despite the QAA's great potential, the 
decade since its introduction has seen the discovery of some 
limitations as a matter of principle. The motivation for this 
work is to understand these limits systematically, clarifying 
where the algorithm will or will not work. This question is of 
broad interest: like in the case of CSA, different failure modes 
point to interesting underlying physics. For instance, critical 
slowing down and the onset of glassiness are two non-trivial 
phenomena which can frustrate CSA. 

Fundamentally, the QAA fails whenever the adiabatic evo- 
lution encounters an exponentially small Hamiltonian gap. 
It is thus tempting to connect the behavior of the adiabatic 
algorithm with various kinds of quantum phase transitions, 
where it is well known that the Hamiltonian gap must close 
in the thermodynamic limit. A simple heuristic suggests that 
first order quantum phase transitions are especially problem- 
atic: the matrix elements of any local Hamiltonian between 
macroscopically distinct states will be exponentially small, 
and hence also the gap of the (barely avoided) crossing. 

Here, we analyze some simple one dimensional first or- 
der transitions and offer both good news and bad news for 
the QAA. The good news is that first-order transitions can be 
accompanied by a finite- size gap which vanishes only alge- 
braically. This is possible because the Hamiltonian gap is not 
a thermodynamic quantity, and is therefore not necessarily en- 
slaved to the details of the transition, in the same way that a 
system exhibiting the breaking of a discrete symmetry can ex- 
hibit gapless excitations on account of frustrating boundary 
conditions. 

The bad news is that we add a failure mode to the known 
limitations of the QAA which we call topological. We con- 
struct a simple classical spin model which has near-degenerate 



ground states which fall into different topological sectors. 
Adding quantum dynamics prefers the sector with exponen- 
tially many ground states, while any degeneracy-lifting inter- 
action favours another containing only 0(1) states. The QAA 
selects the wrong sector in an order-by-disorder mechanism, 
out of which tunneling becomes exponentially slow as the 
quantum dynamics is switched off. 

This is remarkable as our example is a translation invariant 
quasi-one-dimensional Ising model with nearest-neighbour 
interactions only, the ground state of which is readily found 
by inspection, CSA or transfer matrix. This complements rig- 
orous work on the difficulty of finding the precise ground state 
energy of translation invariant ID local Hamiltonians |3J, a 
task which is QMAexp -complete. Of course, the ground state 
of our classical Hamiltonian can be found easily, but its strong 
boundary condition dependence and extreme sensitivity to pa- 
rameters near the quantum first order transition reflect the fea- 
tures that we believe would arise in any QMAexp -complete 
simplification of the Hamiltonians considered in | 3 1. 

Before turning to our results, we briefly review known fail- 
ure modes of the QAA from a physics perspective. Quan- 
tum first order transitions which provably frustrate the QAA 
arise in non-local optimization problems whose energy func- 
tions do not provide "basins of attraction" suitable to local 
exploration in configuration space |4-7|. This reflects the 
inability of local quantum fluctuations to explore non-local 
landscapes effectively - either due to extensive disorder (HO 
or because of golf-course like landscapes with exponentially 
small holes |6, 7|. In models with local energetics, the situa- 
tion is less clear. Thermodynamic calculations within quan- 
tum cavity theory | 8 | suggest random quantum first order 
transitions persist in at least some local models | 9 |, although 
QMC data is inconclusive ifTOll . Recent controversial work 
suggests that Anderson localization may arise in configura- 
tion space when quantum fluctuations are very weak, leading 
to 'perturbative crossings' and exponentially small gaps l lTT- 
14 1 . Heuristic arguments assuming the presence of 'cluster- 
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ing' of pure states in a glassy phase suggest that such cross- 
ings may arise throughout an extended regime of the adiabatic 
evolution ifTSl ; such crossings may even have been observed 
in QMC 1 16 |. In disordered, geometrically local optimization 
problems, Griffiths-like effects may arise in which large local 
regions order before the whole ifTvlfTSll . 

First-order transition with algebraically small gap. Con- 
sider an antiferromagnetic Ising chain in staggered tilted field 
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where J > 0, hi = (— l)*/iisa staggered longitudinal field 
and r is a uniform transverse field. In the thermodynamic 
limit at h = 0, the system exhibits a paramagnetic phase at 
r > J and a Neel ordered phase at F < J where the stag- 
gered magnetization (— exhibits long range order. For 
F < J, as the longitudinal field h is swept down through 
0, the Neel moment rrih = (— l)*(crf )^ jumps from mo+ to 
mo- = — mo+ by a finite amount. This corresponds to a first 
order quantum phase transition where the ground state energy 
density exhibits a first order cusp (as L oo): 



(2) 



Each of the staggered phases exhibits a bulk gap to 
the creation of minority domains. At h = 0, the 
single-wall excitation spectrum may be obtained exactly by 
fermionization, leading to the gapped dispersion e{q) = 
V^J^ + F2 + 2JFcos(7. 

The scaling of the many-body gap at the transition point 
h = at finite size L with periodic boundary conditions 
may be found precisely using fermionization (see Appendix). 
Here, we provide an intuitive perturbative argument near F = 
0+. For L even, the two degenerate Neel ordered ground 
states of H, |1) = lUU - - - U) and |2) = \im • • • It), 
each have energy —LJ. These are separated by a gap 4 J from 
2(2) states with a pair of domain walls. The transverse field 
produces domain walls in pairs and hops individual domain 
walls. In degenerate perturbation theory, the leading order 
at which the transverse field couples within the ground state 
space is L by producing a pair of domain walls and dragging 
them around the system to annihilate. This leads to an expo- 
nentially small gap. 

For L odd, the Neel states do not fit; the lowest energy states 
must have a broken bond somewhere. This can be in any of L 
positions and overall Ising symmetry leads to 2L degenerate 
ground states with energy —J{L — 1) + J: 
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The transverse field acts directly within this state space by 
hopping the frustrated bond : left or right with amplitude — F. 



Thus, the effective Hamiltonian in the ground state space is 
that of a periodic hopping chain of length 2L: 
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The spectrum is Ck = — 2Fcos(/c) with k = |jn and n G 
{0,1,-- - ,2L— 1}. This leads immediately to a unique ground 
state energy of— LJ + 2J — 2Fat/c = with a gap to the 
first excited state of — 2F(cos(^) — 1) ^tt^/L^ which is only 
algebraically small. 

The lesson of this simple calculation is that the Hamiltonian 
gap is not a thermodynamic quantity, in the same way that an 
(in)appropriate choice of boundary conditions can force gap- 
less excitations on a system in which only a discrete symmetry 
is broken. Such non-bulk modes vanish in observables in the 
thermodynamic limit but modify the many-body gap at finite 
size. 

We next present the topological failure mode of the QAA, 
for which we first introduce a quantum dimer model, where 
the topological structure leading to an exponentially small gap 
is most evident. We then transcribe this result to a system 
with an unconstrained Hilbert space - a simple translationally 
invariant Ising ladder. 

Dimer model- Consider a dimer model on a periodic two 
leg ladder of length L. The dimer Hilbert space is spanned 
by hardcore dimer coverings of the ladder. These fall into 
three sectors which are topological in that they are not con- 
nected by any local rearrangement of the dimers. The sectors 
are labelled by a winding number w, the difference between 
the number of dimers on top and bottom rows (on any fixed 
plaquette). On an even length ladder, there are three sectors 
w = ±1,0, while an odd length ladder is topologically triv- 
ial, w = 0. The two 'staggered' sectors, w = ±1 comprise 
only one state each while the = ('columnar') sector has 
2Fl -\- Fl_i ^ states, with Fl the L'th Fibonacci number 
and r = the golden mean. 

The extraordinarily easy search problem we pose to the 
QAA is to find the ground state of the Hamiltonian which 
favours the staggered configurations equally and extensively 
over all of the columnar ones. This is encoded by the local, 
classical Hamiltonian 



^d = ^El')('l+2t/^|=)(= 



(4) 



assigning extensive energy UL to every state in the = 
sector while leaving the two staggered states w = ±1 with 
energy 0. 

It may be surprising that such a simple local Hamiltonian 
can generate a golf-course energy landscape but it is clear that 
such a landscape is hard for the QAA to maneuver. In fact, the 
problem is worse than for a golf course, as we explain next. 

Any off-diagonal term in the dimer Hilbert space involving 
only a finite number of rungs in the ladder leaves the wind- 
ing number w invariant, and hence does not permit transitions 
between the ground states in different winding sectors. Such 
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FIG. 1 . Spectrum of dimer model on even length periodic ladder with 
dimer configurations overlaid. The two staggered w — ±.1 states are 
at ^ = 0, while the columnar sector it; = splits into a band of 
resonating dimer states under resonance F. 



off-diagonal terms induce extensive resonance energies in the 
entropically dominant columnar sector while leaving the en- 
ergy of the lonely staggered states completely unperturbed. In 
other words, not only will switching on a strong quantum per- 
turbation to Hci select the wrong state; switching it off will 
fail to find the correct sector. 

To make this concrete, consider the simplest quantum reso- 
nance term, which flips dimers around a single plaquette: 
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This splits the = sector centred at extensive energy U L 
into a band of width ~ 26LF where 6^0.6 (from exact diag- 
onalisation numerics, which converges rapidly). The ground 
state for w = 0, \RVB), has finite-size energy Ervb = 
UL — BlTL, where Bl converges rapidly to 6 as L ^ oo. 

At a finite resonance field Tc = U /b, \RVB) undergoes a 
strict (unavoided) level crossing with the staggered zero en- 
ergy states (Fig.[T]). In the thermodynamic limit taken through 
even length periodic chains, this crossing corresponds to a first 
order quantum phase transition driven by the resonance: 
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For an odd-length chain, there is only thQw = sector and 
thus no transition at Fc (dreo = —b). This boundary condition 
dependence of the thermodynamic limit reflects the nonlocal- 
ity of the dimer Hilbert space. We address this by transcribing 
the problem into a frustrated Ising ladder in a field. 

Equivalent Ising ladder. The dimer model on the two-leg 
ladder can be turned, via a duality transformation, into a frus- 
trated Ising ladder, the ground states of which map onto the 
dimer states |19t i20il . The simplest Ising model which effec- 
tively reproduces the dimer model physics described above 
turns out to be a two-leg ladder, which features an external 
field and nearest-neighbour interactions only, see Fig. [2] The 
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FIG. 2. The fully-frustrated Ising ladder. Location of Ising spins 
are denoted by filled (empty) circles, subject to fields of strength 
—K {U/2), respectively. Black solid (dashed) lines represent 
(anti-)ferromagnetic interactions of strength K. Dimers live on the 
quasi-dual lattice (fat-dashed lines), with a dimer on the top row de- 
noting an unsatisfied field term. The ground states of the K terms are 
the dimer states, while the U term selects the staggered ones among 
them. 



corresponding classical Hamiltonian is 



Hi = -Y, K,jata] - K + - (7) 

<ij> 



I upper 



i lower 



where Kij = {—)K for solid (dashed) links ij, respectively. 

The ladder is constructed such that (at least) one term of 
order K (spin interaction or field in the top row) per square 
must remain unsatisfied; denoting this term by a dimer placed 
on a fat-dashed line yields the dimer states discussed above. 
Ising configurations with extra dimers are defective; they have 
higher energy in units of K. We observe that non-defective 
staggered dimer configurations correspond to all bottom row 
spins negative. These spins satisfy the U -field, while those in 
the iL' = configurations do not. Thus, Hi encodes the same 
low energy dimer energetics as Hd in the pure dimer model. 

Local quantum dynamics with matrix elements within the 
flippable {w = 0) sector, such as a transverse field F ^ a^, 
lead to the order-by-disorder selection discussed above. This 
favors a resonating state in the = sector over the rigid 
staggered states favored hy U. As the Hilbert space of this 
model is now (local and) larger than for the dimer model, the 
excited state spectrum is richer and the transition more com- 
plicated. The presence of local topological defects at energy 
0{K) breaks the strict conservation of w seen in the dimer 
model, lifting the level crossings. However, for large K/U, 
the first order transition seen in the pure dimer model persists, 
acquiring an exponentially small gap due to the virtual hop- 
ping of defects around the system. 

At large K ^ the pure dimer states are energetically 
well separated from the states with defects. To leading order 
for the pure dimer w = {) sector, the transverse field F on the 
top row acts exactly as Hrvh of Eq. ([5]), while the transverse 
field on the bottom row produces defects. Thus, for infinite K 
and even L, we precisely recover the dimer model and its first 
order transition at Fc = U/b. At large but finite the transi- 
tion is shifted by the resonance energy associated with virtual 
defect states. Thus, to second order, the first order transition 
curve shifts to Fc = U /b ^ /AKb^, in agreement with the 
exact diagonalization results (see Fig. [3]). For odd lengths, the 
first order transition persists although the finite size effects are 
more complicated; nonetheless, the minimum gap arises due 
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FIG. 3. Phase diagram of frustrated two leg ladder. Points indicate 
position of minimum gap in /c = sector as a function of F at fixed 
K, L — 12. Dotted line F — K indicates second order condensation 
predicted to leading order in K/U; dashed line indicates first order 
transition predicted to second order in U/K (see text), (inset) Expo- 
nential scaling of the gap for K in the first order transition regime. 
Dashed lines are best fit exponentials to minimum gap data. From 
top to bottom, K/U = 2.5 - 5.0 in steps of 0.5. 

to the production and dragging of defects around the system 
and it remains exponentially small. 

For K <^ U, the transition from staggered to transla- 
tion invariant phase persists as a second order condensation 
transition. At U K^T, the U field pins the bottom 
row of spins to be negative, leading to an effective Hamil- 
tonian for the top row of an antiferromagnetic Ising chain 
with couplings K in a transverse field F and longitudinal 
field heff = K{1 + (a^J) ^ 2KT'^ /U'^ at leading order. 
This chain indeed exhibits the usual second order transition at 
T ^ K - 0{hlff/K) = K - 0{Kr^/U^) |2F|. In dimer 
language, the liquid state at low K corresponds to a conden- 
sate of defects while that at higher K is 3. resonating dimer 
liquid. 

Exact diagonalisation of the Ising model We back up the 
above assertions with a detailed study of the two leg Ising lad- 
der, for which we describe an interesting phase diagram. We 
exactly diagonalize the zero momentum sector at even lengths 
L from 4tol2(/7 = 1). The points in the phase diagram 
in Fig. |3] are determined by the F corresponding to the mini- 
mum gap at fixed K at size L = 12. Away from the multi- 
critical point where the first order and second order transition 
curves meet, this estimate is indistinguishable from the esti- 
mate made by finite size scaling estimates; near the meeting 
point such estimates are suspect at these small sizes. 

Viewing the two leg ladder Hamiltonian as a optimization 
problem for the QAA, the important feature is the scaling of 
the minimum gap A as a function of L. These are plotted in 
the inset to Fig. |3] for > 2.5 where the transition lies on the 
first order curve. As expected, the exponential fit to A(L) is 
very good, exhibiting a decay constant growing linearly with 
log(iir/F) for K up to 8 (fit not shown). 



Algorithmic implications. It is instructive to consider the 
behavior of the classic "black box" optimization algorithm: 
classical simulated annealing (CSA) with local spin flip dy- 
namics. Similar to the quantum order-by-disorder effect, high 
temperature favors the entropy of the columnar sector, which 
freezes out at low temperatures due to the energetic preference 
for the staggered sector. Unlike the QAA, however, CSA may 
dissipate energy, so that defects in the low temperature regime 
anneal out diffusively with diffusion constants depending on 
temperature but not system size. This leads to 0{L'^) upper 
bounds on the time needed to find the ground state. 

Hence, CSA finds a ground state in the simplest incarna- 
tion. There are ways to fix the QAA as well, but these appear 
to require insight into the nature of the solution and are in 
spirit equivalent to finding the solution by knowing something 
about it, i.e. solving the problem by inspection. 

Conclusions.- The quantum adiabatic algorithm holds 
much promise as a practical tool for solving hard optimiza- 
tion problems, and as such has suffered a barrage of theoreti- 
cal attacks over the last decade. Here, we contribute two sim- 
ple cautionary tales to its analysis: first, that the finite size 
scaling of the Hamiltonian gap exhibited at thermodynamic 
quantum phase transitions must be treated very carefully as 
it may be exponentially sensitive to non-thermodynamic de- 
tails and second, that topology and entropy can overwhelm 
local quantum dynamics, even in translation invariant, local 
qubit systems. Our results also raise the intriguing possibility 
of a further classification of quantum first order transitions by 
non-thermodynamic criteria such as the sensitivity to bound- 
ary conditions of their Hamiltonian gaps. 

Acknowledgments.- We thank E. Farhi, S. Gutman, J. 
Goldstone, D. Gossett and P. Young for discussions. S.L.S. 
was supported by NSF grant number PHY- 1005429. C.R.L. 
acknowledges support from a Lawrence Gollub fellowship 
and the NSF through a grant for ITAMP at Harvard Univer- 
sity. A.S. thanks the CTP at MIT for hospitality and financial 
support. 



[1] E. Farhi, J. Goldstone, S. Gutmann, J. Lapan, A. Lundgren, and 

D. Preda, Science 292, 472 (2001). 
[2] E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser, arXiv 

quant-ph (2000), quant-ph/0001106 
[3] D. Gottesman and S. Irani, FOGS '09. 50th Annual IEEE Sym- 
posium on , 95 (2009). 
[4] E. Farhi, J. Goldstone, D. Gosset, S. Gutmann, and P. Shor, 

Quant. Inf. and Comp. 11, 840 (2011). 
[5] T. Jorg, F. Krzakala, J. Kurchan, and A. C. Maggs, Phys. Rev. 

Lett. 101, 147204 (2008). 
[6] E. Farhi, J. Goldstone, S. Gutmann, and D. Nagaj, Int. J. of 

Quant. Inf. 6, 503 (2008). 
[7] W. van Dam, M. Mosca, and U. Vazirani, in Foundations of 

Computer Science, 2001. Proceedings. 42nd IEEE Symposium 

on (2001) pp. 279-287. 
[8] C. Laumann, A. Scardicchio, and S. L. Sondhi, Phys. Rev. B 

78, 134424 (2008). 



5 



[9] T. Jorg, R Krzakala, G. Semerjian, and F. Zamponi, Phys. Rev. 

Lett. 104, 207206 (2010). 
[10] A. P. Young, S. Knysh, and V. N. Smelyanskiy, Phys. Rev. Lett. 

101, 170503 (2008). 

[11] S. Knysh and V. N. Smelyanskiy, (2010),larXiv: 1005.301 1| 
[12] E. Farhi, J. Goldstone, D. Gosset, S. Gutmann, H. B. Meyer, 

and P Shor, Quant. Inf. and Comp. 11, 181 (2011). 
[13] B. Altshuler, H. Krovi, and J. Roland, PNAS 107, 12446 

(2010). 

[14] M. H. S. Amin and V. Choi, Phys. Rev. A 80, 062326 (2009). 
[15] L. Foini, G. Semerjian, and F. Zamponi, Phys. Rev. Lett. 105, 
167204 (2010). 

[16] A. P. Young, S. Knysh, and V. N. Smelyanskiy, Phys. Rev. Lett. 

104, 020502 (2010). 
[17] B. Reichardt, STOC '04: Proceedings of the thirty-sixth annual 

ACM symposium on Theory of computing (2004). 
[18] D. S. Fisher, Phys. Rev. B 51, 6411 (1995). 
[19] R. Moessner, S. L. Sondhi, and P Chandra, Phys. Rev. Lett. 84, 

4457 (2000). 

[20] R. Moessner and S. L. Sondhi, Phys. Rev. B 63, 224401 (2001). 
[21] A. Ovchinnikov, D. Dmitriev, V. Krivnov, and V. Cheranovskii, 

Phys. Rev. B 68, 214406 (2003). 
[22] A. Kitaev, in Exact Methods in Low -Dimensional Statistical 

Physics and Quantum Computing, edited by J. Jacobsen, S. Ou- 

vry, V. Pasquier, D. Serban, and L. F. Cugliandolo (Oxford 

University Press, Oxford, UK, 2010) Chap. 4. 



Many-body gaps in quantum Ising chain by fermionization 

We consider the quantum Ising chain as in the main text, but 
using a rotated basis of Pauli operators (a^ ^ a^) to connect 
with the convention of |22|: 

i/ = -J^afaJ-r^af (8) 

{ij) i 

We define the IVEajorana Jordan- Wigner transformation: 

j<i 

j<i 

where the 2L operators 7i, 7- are IVEajorana fermion operators 
satisfying the standard anticommutation relations: 

{7i, 7i} = 2<5i,- {7^, 7-} = 25ij {ju 7^} = 
Finally, define the string operators 

Si = l[(-'ljl',)=l[^] (9) 

j<i j<i 

and the total parity operator 

p = sl+1 = n(-nn;) = n (10) 

j 3 



The parity operator P implements the global Ising symmetry 
of the model. 

On a periodic chain of length L, the Hamiltonian may be 
rewritten in the fermionic language 

L-l 

H = -jJ2 5ni5,-+i7i+i - J5L7L71 - r Y.(-'^nj) 

L-l 

j=i 3 

Within each parity sector, the Hamiltonian is quadratic in Ma- 
jorana fermions and corresponds to a bipartite hopping chain 
with a two site unit cell and either periodic (P = — 1) or anti- 
periodic (P = 1) boundary conditions. We may diagonalize 
such a problem by Fourier transform to find (positive energy 
modes): 

e{k) = Vj2+r2-2jrcos(A:) (11) 

where we take /c ^nforP -1 and/c ^(n + 1/2) for 
P = l,n = 0- -- P — 1. The vacuum state for the quadratic 
problem H± with P = ±1 is ±) and has energy 

±) = - ^ e{k) \n, ±) = El \n, ±) (12) 

k 

where k runs over the appropriate momenta for the boundary 
condition. 

In Eq. ([12]) we have calculated the energy of ±) within 
the effective Hamiltonian H± which fixes the value P = ±1, 
ignoring the consistency condition that actually have 

parity P. For a general quadratic Majorana Hamiltonian, 

H = ^Jihijjj (13) 

where hij is a real antisymmetric coefficient matrix and z, j 
run over 1, 1', 2, 2' • • • , the lowest energy state \Q) has parity 

P\Q)=sgn{Ff{h))\Q). (14) 

If P ±) = ± ±), the vacuum energy corresponds to 
the ground state energy in that parity sector. If not, the ground 
state in the sector will correspond to cj ±) where cj cre- 
ates the lowest energy excitation of the quadratic Hamiltonian 
as this is the lowest state with the correct parity. 

To compute the parity of the vacua for the magneti- 
cally ordered phase, we exploit the adiabatic invariance of 
sgn(Pf(/i)). Thus, we may work at the trivial points in the 
phase diagram where either J or F is 0. First, we observe that 
for J = 0, F = 1 in the paramagnetic phase, the Hamiltonian 
reduces to 

H = iY.^iii (15) 

whose ground state clearly sets —ijj'Jj = 1 for all j. Thus 
the parity of the vacuum is 1 for both boundary conditions 



6 



and the system has a unique ground state in the P = 1 sector 
with a gap to single particle excitations in the P = —1 sector 
as expected. 

Now, we consider the ferromagnetic J = 1 , F = Hamil- 
tonian: 

H = iY,l'jlj+i-P^l'Lli (16) 

For P = —1, hpM for this Hamiltonian is precisely a 
translation by 1 step in the length 2L chain of Majoranas. 
Thus, hpM = RhpuR^ where R translates every site by 
1. Since Pf (mpMi^^) = det (i^) Pf (/ipm), and det(i?) = 
(-1)2^-1 = -1, we have that P|n,-) = Simi- 
larly, flipping P to +1 multiplies one row and column of hpM 
by -1, flipping sgn(Pf(/iFM)). Thus, P \VL, ±) = ± ±) in 
the ferromagnetic phase and therefore the ground state of each 
sector is the vacuum state for the effective quadratic Hamilto- 
nian in that sector. Thus, for the ferromagnet, we have the 
bottom two states have energy: 

n 
n 

Since e{k) is analytic and lives on the circle k G (0, 27r), the 



Fourier series coefficients decay faster than any power law 
in m. Using this, it is straightforward to show that the — 
E^_ is smaller than any power law in L. 

Finally, we turn to antiferromagnetic J = — 1 , F = 0. In 
this case, 

hAFM = —hpM (17) 

and thus Pf(/iAFM) = (— 1)^ Pf(^FM)- For L even, this 
means that there are two ground states identical to those found 
in the ferromagnetic case with an exponentially small separa- 
tion in L. For L odd, the parity of the vacuum in both sectors 
is wrong: 

P\n,±) = T\^.±) (18) 

and thus the two lowest energy states correspond to populating 
the lowest energy excitation in each sector: 

^+ = ^^^26(7r+^) (19) 
E_=E^_^ 2e(7r) (20) 

Expanding the single particle dispersion e at its quadratic min- 
imum at A: = TT, this leads to a power law small gap 

^+-^_^2F^. (21) 



